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INTRODUCTION 


■x 

Steady-state  vibration  response  analyses  are  commonly  used  in  evaluating  the 
dynamic  behavior  of  structures  subjected  to  cyclic  external  forces.  An  undocumented 
vibration  response  analysis  based  on  modal  superposition  was  developed  about  10  years 
ago.  The  analysis  calculates  the  acceleration  response  at  any  selected  point  on  a 
structure  for  specified  vibratory  loading;  the  program  provides  tabular  output  of 
modal  accelerations  and  acceleration  response  with  optional  paper  plots. 

In  using  the  program,  it  became  obvious  that  several  improvements  were  desir¬ 
able;  thus,  an  effort  was  undertaken  to  expand  the  program's  capability,  to  decrease 
c'  "'turnaround  time,  "and  to  make  the  program  easier  to  use.  The  mathematical  equations 
used  in  the  original  program  were  retained,  but  a  new  program  logic  that  would  use 
less  memory  storage  was  developed.  The  new  program  operates  interactively  for  rapid 
evaluation  and  plotting  of  structural  response  to  vibratory  loads.  Options  are  also 
provided  to  control  both  the  printed  and  the  interactive  graphic  output. 

t  1  li  ^  ^  \  *  *  rx.f 'c  r  i 

The  purpose  of  this  report  is  to  document  the  improved  version  of  thrs~  computer 
program.  The  theoretical  background,  program  description,  and  application  are  pre¬ 
sented  along  with  user  instructions  and  a  sample  interactive  computer  session^^  The 
program  can  use  results  from  any  vibration  analysis  that  has  output  in  the  form  of 
generalized  masses,  eigenvalues,  and  eigenvectors.  For  instance,  the  program  can  be 
used  as  a  postprocessor  with  finite-element  codes  such  as  NASTRAN®  (ref.  1),  or  it 
can  use  measured  modal  data  (as  discussed  in  ref.  2) .  This  computer  program,  named 
VIBRA  (vibration  response  analysis),  is  available  through  COSMIC1  as  LAR-13291. 


SYMBOLS 

A  modal  acceleration,  g/lb 

[D(w)J  damping  matrix,  lb-sec/in. 

d(u>)  generalized  modal  damping,  lb-sec/in. 

{F (  w)  }  applied  force  vector,  lb 

g  acceleration  of  gravity  (G  in  computer-generated  tables  and  plots),  in/sec 

gn  structural  damping  coefficient  (eq.  (12)) 

i  ; 

i  imaginary  operator 


j  response  coordinate 

(K)  stiffness  matrix,  lb/in. 

k  force  coordinate 


'  r 

r  v  .  V  r 


COSMIC,  112  Barrow  Hall,  Univ.  of  Georgia,  Athens,  GA  30602. 


M 

[M] 

N 


2 

generalized  mass,  lb-sec  /in. 

.2, 


mass  matrix,  lb-sec  /in. 
number  of  modes 
n  mode  number  index 

{q ( o>)  }  displacement  response  vector,  in. 

{q(cj)}  acceleration  response  vector,  g  units 
[Y(<d)  ]  displacement  mobility  matrix,  in/lb 
[Y(  (j)  ]  acceleration  mobility  matrix,  g/lb 
tZ(w)]  displacement  impedance  matrix,  lb/in. 

[Z(u)>]  acceleration  impedance  matrix,  lb/g 

C  critical  damping  ratio 

<J>  normalized  mode 

<|»  orthonormal  mode,  in1  /2/lb1  /2-sec 

w  forcing  frequency 

u>n  natural  frequency  of  nth  mode,  Hz 

Superscripts: 

I  imaginary  component  of  conplex  variable 

R  real  component  of  complex  variable 

T  matrix  transpose 

Dots  over  a  symbol  define  parameters  in  terms  of  acceleration. 


METHOD  OF  ANALYSIS 

The  matrix  equations  of  motion  describing  the  dynamic  behavior  of  linear,  pro¬ 
portionally  damped  structures  can  be  written  in  the  frequency  domain  (ref.  2)  by 
assuming  simple  harmonic  motion: 


(  —  tot  M  J  +  iw[D{(0)]  +  [  K]  )  {q  (  oj)  }  =  {F  (  w)  } 


(1  ) 


2 


V  S  *».  T — 


(u>/io  )  d  (w) 
n  n 


2  2 
+  d  (uj) 
n 


where  co^  is  the  natural  frequency  of  the  nth  normal  mode  and  is  the  (jk)th 

modal  acceleration  of  the  nth  normal  mode  defined  by 


in  kn 


The  orthogonal  modal  vector  {$}  and  mass  matrix  [M]  are  combined  to  calculate 
the  generalized  mass  Mn  in  equation  (8)  by 


{<)>}  =  M 

n  n  n 


If  the  modal  vector  is  scaled  such  that 


U)  [  M  ]  { <(» }  =  1.0 

n  n 


then  the  modal  vector  is  the  orthonormal  modal  vector  (4>)  ,  and 

n 

A  =  4>.  <1\ 
jkn  jn  kn 

Reference  2  includes  a  complete  derivation  of  the  modal  superposition  equations 
described  in  equations  (6)  and  (7).  The  modal  damping  d  (w)  can  be  defined  as 
structural  or  viscous.  For  structural  damping  n 


d  (to)  =  g 
n  n 


and  for  viscous  damping 


d  (oj)  =  2C  —  (12 

n  a) 

n 

where  C  is  the  critical  damping  ratio. 

In  principle,  the  number  of  modes  N  is  infinite;  in  practice,  only  a  finite 
number  of  modes  are  necessary  over  a  specified  frequency  range.  For  unconstrained 
structures,  N  includes  both  rigid-body  and  elastic  modes.  These  modal  parameters 
can  be  derived  from  either  test  or  analysis,  as  appropriate. 
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where  [M J  and  [ K ]  are  real  and  symmetrical  mass  and  stiffness  matrices,  and  the 
symmetrical  matrix  [D(oj)]  represents  a  general  form  of  damping.  The  terms  on  the 
left  side  of  equation  (1)  define  the  displacement  impedance  matrix  [ Z ( co ) ]  such  that 


tz (to )  Hq(co)}  =  {F(w)} 


Similarly,  the  acceleration  impedance  [z ( to) 1  is  defined  as  the  matrix  coefficients 
of  accelerations,  and  equation  (2)  becomes 


(Z(u)){q(u)}  =  { F(w» 


The  acclerations  are  obtained  from  equation  (3)  by 


{q(w )}  =  [Y(io )  ]  {  F(w)} 


where 


[Y(u))J  =  tZ(to)] 


The  variables  q  and  Y  are  not  differentiated  with  respect  to  time  but  are 
intended  to  define  responses  in  terms  of  acceleration.  The  variables  in  equation  (4) 
are,  in  general,  complex  valued  and  frequency  dependent.  Ihe  matrix  [Y((j)]  defined 
as  acceleration  mobility,  is  a  transfer  function  which  relates  input  excitations 
{F  ( w)  }  to  output  accelerations  {q(u))}.  Acceleration  responses  are  used  in  this 
formulation  because  test  measurements  are  generally  obtained  from  acceleration  trans¬ 
ducers.  However,  the  acceleration  mobility  Y  can  be  related  to  the  displacement 
mobility  Y  in  the  frequency  domain  by 


Y(u>)  =  -w  Y(u>) 


Physically,  the  (jk)th  element  in  the  acceleration  mobility  matrix  Y  defines 

3k 

the  response  at  j  due  to  a  force  or  moment  at  k.  The  acceleration  mobility  is 

.  ••  ..I 

conposed  of  real  and  imaginary  components,  where  Y  =  Y  +  Y.  ,  and  can  be 

expressed  in  terms  of  modal  parameters:  3  3  3k 


Y*  = 
3k 


(w/U)  )  1  -  (oj/co  ) 

_ n  l  _ n  j 

1  -  <w/u>  )2  2  +  d  (0,)2 
n  n 
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Appendix  C  contains  information  on  computer  memory  requirements,  auxiliary  storage 
files,  and  graphics  subroutines. 


APPLICATIONS 

Applications  are  presented  for  a  constrained  viscously  damped  system  (case  1 ) 
and  an  unconstrained  structurally  damped  helicopter  tail  boom  (case  2).  These  appli¬ 
cations  were  selected  to  demonstrate  the  versatility  of  the  VIBRA  program. 

Case  1  -  Constrained  Viscously  Damped  System 

Case  1  demonstrates  two  methods  of  applying  forces  to  a  system.  The  modal  data, 
given  in  table  I,  are  for  the  constrained  system  with  no  damping.  However,  the  damp¬ 
ing  terms  may  be  included  when  calculating  response  because  of  the  modal  superposi¬ 
tion  principle. 

Case  1  A,  as  shown  in  figure  2(a),  is  the  constrained  system  with  a  unit  force 
applied  at  node  1  with  response  to  be  calculated  at  node  1  (the  driving  point)  and  at 
node  5  (a  transfer  point),  assuming  2.5  percent  viscous  damping  for  all  modes. 

Table  II  shows  the  VIBRA  input  data  file  7.  The  problem  title  and  the  control  param¬ 
eters  given  in  line  2  of  table  II  are  described  in  table  III.  The  force  has  a  fre¬ 
quency  ranging  from  0  to  50  Hz  with  an  increment  of  1  Hz.  Acceleration  responses  for 
nodes  1  and  5  are  given  in  tables  IV  and  V,  while  the  corresponding  plots  are  shown 
in  figures  3  and  4. 

Case  IB,  as  shown  in  figure  2(b),  has  multiple  forces  applied  at  nodes  1,  2, 
and  3,  assuming  variable  viscous  damping.  Tables  VI  and  VII  describe  the  force  vec¬ 
tor  sets  and  damping  data,  respectively.  The  force  vectors  are  composed  of  real  and 

imaginary  parts  with  different  magnitudes  and  directions.  Both  vector  sets  have  a 

20-Hz  forcing  frequency.  Because  the  responses  are  calculated  for  only  one  fre¬ 
quency,  no  plots  can  be  generated.  The  input  data  file  for  case  IB  is  given  in 
table  VIII,  with  the  problem  title  and  control  parameters  description  shown  in 

table  IX.  Responses  are  calculated  at  all  locations  for  both  vector  sets.  Tables  X 

and  XI  show  the  modal  acceleration  matrix,  acceleration  mobility,  and  acceleration 
response  for  nodes  1  and  2.  Table  XII  gives  the  acceleration  response  for  nodes  3, 

4,  and  5. 


Case  2  -  Unconstrained  Structurally  Damped  Helicopter  Tail  Boom 

Undamped  modal  data  for  the  helicopter  tail  boom  shown  in  figure  5  were  obtained 
from  a  finite-element  analysis  and  reformatted  for  input  to  VIBRA.  The  assumed  damp¬ 
ing  terms  are  included  in  the  input  file.  As  figure  5  shows,  a  unit  force  is  applied 
at  the  end  of  the  tail  boom  (node  location  193)  and  the  response  is  calculated  near 
the  midsection  (node  109).  Table  XIII  shows  the  input  data  file  7.  Only  modal  data 
at  the  load  and  response  node  locations  are  required.  Data  are  input  for  a  total  of 
12  orthonormal  modes  of  which  the  first  6  are  rigid-body  modes.  The  forcing  fre¬ 
quency  ranges  from  75  to  300  Hz  with  a  frequency  increment  of  1  Hz.  Since  the  modes 
are  orthonormal,  the  generalized  mass  has  a  value  of  1.0  for  all  modes.  The  assumed 
structural  damping  coefficients  are  input  on  line  10  of  table  XIII,  with  values  rang¬ 
ing  from  1.5  to  1 .6  percent.  The  problem  title  and  control  parameter  description  are 
given  in  table  XIV. 


PROGRAM  ORGANIZATION 


Program  Description 

The  VIBRA  program  (vibration  response  analysis)  uses  equations  (6)  and  (7)  to 
develop  the  acceleration  mobilities  and  equations  (8)  and  Ml)  to  form  the  modal 
accelerations.  The  structural  responses  are  determined  for  known  force  excitations 
using  equation  (4).  Therefore  the  structure  to  be  analyzed  must  be  described  by  a 
consistent  set  of  eigenvectors,  eigenvalues,  and  generalized  mass  data.  These  modal 
data  may  be  obtained  from  any  natural  vibration  analysis  or  from  ground  vibration 
measurements  of  the  structure.  Only  the  modal  data  at  the  selected  response  and 
applied  vibratory  load  points  are  required  to  calculate  the  forced  response  of  the 
structure. 


Program  Flow 

A  flow  chart  of  the  VIBRA  program  is  shown  in  figure  1  with  prepared  input  data 
read  from  file  7.  As  shown  in  the  figure,  the  program  reads  in  a  title  card,  a  con¬ 
trol  card,  natural  frequencies  w  ,  damping  coefficients  Q,  generalized  mass  M^, 
and  one  or  more  force  vector  sets11  {F(w)}  followed  by  the  modal  coordinates 
{4>n=i  at  the  respective  load  point  k.  The  complex  valued  force  vector  {F(u))} 

defines  the  magnitude  and  direction,  while  the  modal  coordinates  {<(>  }  at  the 

applied  load  point  k  define  the  location  of  the  force  vector.  Next,  one  set  of 
modal  coordinates  {$n=1  N  }j  at  the  response  point  j  is  read  from  input  file  7. 

The  program  calculates  one  response  point  at  a  time.  The  modal  accelerations,  accel¬ 
eration  mobilities,  and  acceleration  responses  are  calculated  for  this  response 
point.  Then,  print  options  are  available  to  output  these  results  either  on  the  ter¬ 
minal  screen  or  output  file  8,  or  both.  Next,  VIBRA  produces  a  plot  which  displays 
the  acceleration  response  q  versus  forcing  frequency  u>  in  four  component  plots 
consisting  of  (1  )  real  acceleration  versus  go,  (2)  imaginary  acceleration  versus  go, 
(3)  amplitude  versus  go,  and  (4)  phase  angle  versus  co.  After  reviewing  the  plot, 
the  user  has  the  option  to  plot  any  of  the  conponent  plots  separately  on  the  screen 
or  continue  execution.  At  this  point  the  user  specifies  a  value  for  the  NAMELIST 
VPLOT  parameter  NPLOT  to  indicate  the  next  desired  operation. 

Figure  1  shows  the  program  logic  determined  by  the  value  specified  for  NPLOT.  A 
positive  value  of  1,  2,  3,  4,  or  5  for  NPLOT  causes  the  program  to  display  the  plot 
and  then  continue  with  the  next  desired  program  operation.  For  a  negative  value,  -2, 
-3,  -4,  or  -5,  no  plot  is  created  but  the  program  continues  with  the  next  desired 
operation.  A  value  of  NPLOT  =  1  instructs  the  program  to  read  another  NAMELIST 
VPLOT.  A  value  of  NPLOT  =  +2  or  -2  causes  the  program  to  apply  the  next  force  set, 
if  there  is  more  than  one  input,  to  the  present  response  case.  A  value  of 
NPLOT  =  +3  or  -3  instructs  the  next  response  group  to  be  read  from  input  file  7.  If 
the  previous  response  case  was  calculated  using  other  than  the  first  force  vector 
set,  the  first  force  vector  set  will  be  used  again  when  the  next  response  case  is 
read.  A  value  of  NPLOT  =  +4  or  -4  causes  the  program  to  read  a  title  card  for  a 
new  problem.  To  exit  the  program,  a  value  of  NPLOT  =  +5  or  -5  must  be  used. 

Appendix  A  contains  detailed  information  on  the  data  and  format  requirements  for 
the  prepared  input  file  7.  The  type  of  interactive  input  required  by  VIBRA  and  a 
description  of  the  namelist  parameters  are  provided  also.  To  demonstrate  the  flow 
and  graphics  capability  of  the  program,  an  example  problem  is  solved  in  an  interac¬ 
tive  session  provided  in  appendix  B.  Also  included  in  appendix  B  is  a  brief  descrip¬ 
tion  of  the  example  problem's  modal  data  and  hew  the  data  are  arranged  in  file  7. 
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A  plot  of  node  109  response  is  shown  in  figure  6;  the  real  acceleration  cotrpo- 
nent  is  examined  further  in  figures  7  and  8  by  using  the  zoom  feature.  The  zoom 
feature,  described  in  appendix  B,  allcws  the  user  to  define  a  portion  of  a  plot  to  be 
expanded.  Figure  8  shows  a  flattened  peak  around  212  Hz.  This  effect  is  due  to  not 
calculating  response  exactly  at  the  natural  frequency.  In  view  of  these  results, 
care  must  be  taken  when  determining  the  frequency  increment.  The  imaginary  component 
of  the  acceleration  response  is  shown  in  figure  9,  and  the  amplitude  response  is 
shown  in  figure  10. 


CONCLUDING  REMARKS 

An  interactive  computer  program  (VIBRA)  for  calculating  steady-state  frequency 
response  of  linear,  proportionally  damped  structures  has  been  described.  The  program 
uses  a  modal  superposition  approach  to  calculate  the  structural  response.  The 
response  is  calculated  as  a  function  of  frequency  and  in  terms  of  acceleration.  The 
program  has  interactive  graphics  capabilities  to  display  the  acceleration  response 
versus  forcing  frequency.  The  interactive  capability  enables  rapid  evaluation  of  the 
structural  response.  Sample  problems  consisting  of  several  simple  systems  and  a 
complex  system  demonstrate  program  versatility. 


Langley  Research  Center 

National  Aeronautics  and  Space  Administration 
Hampton,  VA  23665 
May  29,  1984 


APPENDIX  A 


INPUT  DATA  DESCRIPTION 

VIBRA  Input  File  7  Data  and  Format  Requirements 

The  input  consists  of  prepared  data  on  file  7  and  interactive  input  from  a  key¬ 
board.  VIBRA  reads  the  input  data  file  7  using  list-directed  read  statements. 
Therefore,  all  numerical  data  can  be  entered  on  a  card  image  in  any  column.  The  only 
requirement  is  that  the  data  must  be  separated  by  blanks  or  commas.  A  blank  space 
cannot  be  used  for  a  zero.  As  shown  in  tables  II,  VIII,  and  XIII,  a  typical  input 
file  consists  of 

1 .  A  card  containing  the  problem  title 

2.  A  control  card  specifying  parameters 

3.  Natural  frequency  group 

4.  Generalized  mass  group 

5.  Damping  group 

6.  Force  vector  group(s) 

7.  Force  coordinate  group(s) 

8.  Response  coordinate  group(s) 

A  data  group  is  defined  by  a  single  title  card  with  information  identifying  ie 
group,  followed  by  one  or  more  cards  containing  the  data  values.  The  data  value  are 
entered  in  free-field  format.  Except  for  the  force  vector  group,  where  only  a  l- 
mum  of  six  groups  may  be  input,  there  is  no  limit  on  the  number  of  force  and  re:  _.r 
coordinate  groups  that  can  be  input  for  one  problem.  Also,  there  is  no  limit  on  ^ 
number  of  problems  that  can  be  on  the  same  input  file.  The  description  and  format 
for  the  data  file  are  given  in  the  following  eight  sections. 

Title  card.-  A  single  card  contains  the  problem  title. 


FORTRAN 
variable  name 


Description 


TITLE  Heading  information  for  printout 

Control  card.-  Each  of  the  10  parameters  listed  on  the  control  card  must  be 
specified  and  entered  in  the  following  order: 

FORTRAN 

variable  name  Description 

NZ  Number  of  rigid  body  modes 

NM  Total  number  of  modes  (rigid  +  elastic) 

MNORM  Type  of  mode: 

1  orthonormal  modes  (one  term  input) 

0  normal  modes  (NM  terms  input) 
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FORTRAN 

variable  names  Descripti an 

Type  of  modal  damping: 

1  viscous  damping 

0  structural  damping 

Variable  damping  control: 

1  damping  is  the  same  for  all  modes  (one  term  input) 

0  damping  varies  for  each  mode  (NM  terms  are  input) 

Number  of  force  vector  groups  (maximum  of  6) 

Number  of  forces  applied  to  structure  (if  NF  >  1 ,  then  only 
frequency  point  will  be  calculated  using  the  WSTRT  value) 

Starting  frequency  point,  Hz 

Ending  frequency  point,  Hz 

Frequency  step 

Natural  frequency  group.-  The  natural  frequency  group  is 

FORTRAN 

Card  variable  name  Description 

1  WLABEL  Group  title  card 

2+  WN,  Hz  Natural  frequency  for  nth  mode 

Note:  There  should  be  NM  terms  input. 

Generalized  mass  group.-  The  generalized  mass  group  is 


MDAMP 

NZETA 

NFS 

NF 

WSTRT 

WSTOP 

DELW 


Card 

1 

24- 


FORTRAN 
variable  name 

MLABEL 

M,  lb-sec2/in. 


Description 

Group  title  card 

Generalized  mass  for  nth  mode 


Note:  If  MNORM  =  1,  one  mass  term  is  input  (orthonormal  modes);  If  MNORM 

NM  terms  are  input. 
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Damping  group.-  The  damping  group  is 
FORTRAN 

Card  variable  name  Description 

1  GLABEL  Group  title  card 

2+  G  Modal  damping  coefficient 

(for  viscous  damping,  enter 
for  structural  damping  enter  2 £) 

Note:  If  NZETA  =  1,  one  damping  term  is  input  (damping  same  for  all  modes);  if 

NZETA  =  0,  nm  terms  are  input. 

Force  vector  group.-  The  force  vector  group  is 
FORTRAN 

Card  variable  name  Description 

1  FLABEL  Group  title  card 

2+  F,  lb  Complex  force  vector,  (F^F1) 

Note:  There  should  be  NFS  force  vector  groups  and  NF  forces  for  each  vector  group. 

The  order  in  which  each  force  is  entered  also  determines  the  order  that  the  force 
coordinate  groups  are  input.  The  vector  group  forms  the  magnitudes  and  directions 
of  the  forces  applied  to  the  structure.  A  maximum  of  six  force  vector  groups  may 
be  input. 

Force  coordinate  group.-  The  force  coordinate  group  is 
FORTRAN 

Card  variable  name  Description 

1  PLABEL  Group  title  card 

2+  PHI  Force  modal  coordinate  for  nth  mode 

Note:  There  should  be  NF  force  coordinate  groups  and  NM  terms  for  each  group.  The 

force  coordinate  groups  must  follow  the  same  input  order  as  the  force  vectors.  The 
force  coordinate  defines  the  location  where  the  vectors  are  applied.  When  more 
than  one  force  coordinate  group  is  entered,  only  one  frequency  point  is  calculated 
for  the  response. 


» 


» 


» 


» 
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Response  coordinate  group.-  The  response  coordinate  group  is 


FORTRAN 

Card  variable  name  Description 

1  RLABEL  Group  title  card 

2+  PHIBAR  Response  modal  coordinate  for  nth  mode 

Note:  There  should  be  NM  terms  input  for  each  group.  There  can  be  as  many  groups  as 
desired,  the  only  limit  being  the  total  number  of  degrees  of  freedom  of  the  struc¬ 
ture.  To  calculate  the  response  at  the  driving  point,  input  the  particular  force 
coordinate  group  as  a  response  group  also. 


VIBRA  Interactive  Input 

Input  from  the  interactive  terminal  keyboard  is  prompted  by  the  program.  Key¬ 
board  input  includes  pressing  the  carriage  return  key  to  advance  to  the  next  program 
section,  answering  questions  with  the  character  Y  or  N  for  yes  or  no,  specifying 
print  options  for  the  output  of  results,  and  entering  NAMELIST  VPLOT  data.  In  speci¬ 
fying  print  options,  the  program  provides  a  screen  menu  with  input  instructions. 
NAMELIST  VPLOT  is  used  to  specify  program  flow  and  plot  options.  In  the  first  part 
of  this  section,  the  NAMELIST  parameters  and  their  default  values  are  described;  in 
the  second  part,  NAMELIST  rules  are  given. 

Contents  of  NAMELIST  VPLOT.-  Descriptions  of  the  NAMELIST  parameters  and  their 
default  values  are  as  follows: 

QREAL 

1  Plot  real  acceleration  versus  forcing  frequency 
0  No  plot  (default) 

QIMAG 

1  Plot  imaginary  acceleration  versus  forcing  frequency 
0  No  plot  (default) 

AMPL 

1  Plot  acceleration  amplitude  versus  forcing  frequency 
0  No  plot  (default) 

PHASE 

1  Plot  phase  angle  versus  forcing  frequency 
0  No  plot  (default) 
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U SCALE 

1  Input  minimum  and  maximum  values  for  X-  and  Y-axes 
0  Automatic  program  scaling  of  plot  axes  (default) 

If  U SCALE  =  1,  the  following  parameters  must  be  selecteu  and  specified  to 
change  the  default  values: 

XMIN  is  minimum  value  for  X-axis  (forcing  frequency  value);  default,  0.0 

XMAX  is  maximum  value  for  X-axis  (forcing  frequency  value);  default,  10.0 

YMIN  is  minimum  value  for  Y-axis  (response  value);  default,  0.0 

YMAX  is  maximum  value  for  Y-axis  (response  value);  default,  10.0 

The  following  parameters  may  be  specified: 

JLINE  Selects  number  of  data  points  to  be  plotted: 
n  plot  every  nth  point 
1  (default) 

XDIV  Number  of  minor  tick  mark  divisions  for  X-axis  (default,  10.0) 

YD IV  Number  of  minor  tick  mark  divisions  for  Y-axis  (default,  10.0) 

LTYPE  Type  of  line  for  plotting: 

-1  no  line 

0  solid  (default) 

1  long  dash 

2  dash 

3  dot  dash 

4  dot 

ISYM  Symbol  type: 

0  no  symbol  (default) 

1  circle 

3  triangle 

4  square 

5  star 

6  diamond 

7  plus  sign 

ISIZE  Multiplicative  factor  for  symbol  size: 

1  small  symbols  (default) 

2  medium  symbols 

3  large  symbols 
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NPLOT  Control  loop  parameter: 

1  plot  response  then  read  another  NAMELIST  VPLOT 

2  plot  response  then  apply  next  force  vector  set 

3  plot  response  then  read  another  response  group 

4  plot  response  then  read  another  problem 

5  plot  response  then  exit  the  program  (default) 

-2  no  plot  but  apply  next  force  vector  set 

-3  no  plot  but  read  another  response  group 

-4  no  plot  but  read  another  problem 

-5  no  plot  and  exit  the  program 

Note:  Figure  1  shows  how  the  NPLOT  parameter  controls  looping  within  the  program. 

Also,  all  previous  responses  must  have  been  calculated  to  position  the  input  file 
for  reading  another  problem  title. 

NAMELIST  rules .-  The  following  NAMELIST  rules  apply  to  most  FORTRAN  systems: 

1.  A  NAMELIST  begins  in  column  2  with  a  dollar  sign,  e.g., 

12345678910  .  80 

$VPLOT 

2.  NAMELIST  parameters  are  separated  by  commas. 

3.  NAMELISTS  end  with  a  blank  and  then  a  dollar  sign,  no  final  comma,  e.g., 

$VPLOT  QREAL=1 ,NPL0T=+2  $ 

4.  If  more  than  one  line  is  used  to  enter  a  namelist,  then  that  line  must  end 
with  a  comma,  e.g., 

$ VPLOT  QREAL=1 ,NPLOT=2 , USCALE=1 , 

XMIN=0. 0, XMAX=45. 0, YMIN=0. 0, YMAX=5. 0  $ 


Blank  Common  Length 

VIBRA  uses  the  control  card  parameters  to  determine  the  location  of  program 
variables  within  the  blank  common  array  ZZZ,  which  has  a  set  length  of  5000  decimal 
memory  locations.  The  number  of  required  memory  locations  for  a  problem  can  be 
determined  from 

BCL  =  5* ( 1 +NM )  +  NF* (NM+2*NFS+4 )  +  (5*IFN)  +  6 

where 

IFN  =  (WSTOP  -  WSTRT ) /DELW  +  1 
which  is  the  number  of  forcing  frequency  step  points. 
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SAMPLE  INTERACTIVE  SESSION 

In  the  sample  interactive  session,  responses  are  calculated  at  mass  1  and  mass  2 
of  the  unconstrained  spring  mass  system  shown  in  figure  B1  for  a  unit  force  applied 
at  mass  3.  The  force  has  a  frequency  ranging  from  0  Hz  to  20  Hz  in  0.5-Hz  incre¬ 
ments.  Modal  data  for  the  system,  obtained  from  a  simple  vibration  analysis  program, 
are  shown  in  table  BI.  The  system  has  one  rigid-body  and  two  elastic  normal  modes. 

The  input  data  file  7  for  this  program  is  given  in  table  BII.  The  natural  fre¬ 
quencies  and  generalized  masses  from  table  BI  are  input  in  modal  order  as  shown  in 
table  BII,  lines  4  and  6,  respectively.  The  complex  force,  table  BII,  line  10,  has  a 
real  component  term  only.  In  order  to  apply  the  force  at  mass  3,  the  modal  ampli¬ 
tudes  at  mass  3  are  entered  as  the  force  coordinates,  as  shown  in  line  12  of 
table  BII.  Similarly,  to  evaluate  the  responses  at  mass  1  and  mass  2,  modal  ampli¬ 
tudes  at  those  points  are  entered  in  line  14  and  line  16,  respectively.  If  the 
response  at  mass  3  was  desired,  the  corresponding  modal  amplitudes  would  have  to  be 
reentered  as  response  coordinates  also. 

After  the  input  file  7  is  prepared,  the  user  is  ready  to  execute  VI3RA  on  an 
interactive  graphics  terminal.  This  section  describes  a  typical  interactive  session 
using  the  input  file  7  for  the  spring  mass  system,  table  BII.  For  each  new  problem 
the  first  screen  page  contains  a  title  and  descriptive  table  of  control  parameters, 
as  shown  in  table  Bill.  VIBRA  pauses  after  each  screen  page  to  allow  for  viewing  and 
obtaining  a  hard  copy.  When  the  terminal  bell  rings,  the  program  temporarily  stops 
execution  until  the  return  key  is  pressed.  Pressing  the  return  key  erases  the  screen 
and  causes  VIBRA  to  continue  to  the  next  program  operation. 

After  VIBRA  calculates  a  response  case,  a  print  menu  is  displayed.  As  shown  in 
table  BIV,  option  2  is  selected  for  the  first  three  questions.  This  causes  the  modal 
accelerations,  acceleration  mobilities,  and  response  tables  to  be  output  to  auxiliary 
file  8.  The  next  two  questions  control  print  output.  In  this  example,  50  lines  of 
output  per  page  for  file  8,  using  a  print  step  of  1,  were  selected.  A  print  step  of 
1  causes  every  frequency  point  to  be  printed.  Upon  completion  of  the  print  menu, 
VIBRA  pauses  until  the  carriage  return  is  pressed. 

Next,  a  preview  plot  of  the  response  components  is  given  automatically,  as  shown 
in  figure  B2.  VIBRA  pauses  upon  completion  of  the  plot  until  the  carriage  return  key 
is  pressed,  which  causes  VIBRA  to  prompt  for  NAMELIST  VPLOT  input.  As  table  BV 
shows,  the  real  acceleration  component  plot  is  selected  by  setting  the  QREAL  param¬ 
eter  equal  to  1 .  The  NPLOT  parameter  is  defined  to  be  a  positive  3.  This  means  the 
selected  component  will  be  plotted  and  then  another  response  will  be  calculated  using 
the  same  force  vector  set  as  the  previous  group.  After  the  grid  question  in  table  BV 
'  is  answered,  VIBRA  immediately  erases  the  screen,  and  the  response  is  plotted  in  fig¬ 

ure  B3  with  a  question  at  the  bottom  of  the  screen  to  zoom  in  on  the  graph.  The  zoom 
option  enables  the  user  to  examine  any  portion  of  a  response  plot  on  an  expanded 
scale.  The  user  defines  the  scale  by  entering  two  position  coordinates,  upper  left 
and  lower  right.  Figure  B3  shows  that  the  zoom  option  is  selected.  Cursor  cross 
i  hairs  appear  on  the  screen  and  are  moved  by  using  the  terminal  thumbwheels.  The  two 

|  position  coordinates  are  entered  one  at  a  time.  The  coordinates  are  input  by  posi- 

'  tioning  the  cross  hairs  at  the  desired  point,  pressing  any  key  (followed  by  a  car¬ 

riage  return),  and  waiting  for  the  terminal  bell  to  ring.  The  cross  hairs  disappear 
momentarily  until  the  coordinate  is  received;  then  the  terminal  bell  rings.  After 
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the  coordinates  are  entered,  VIBRA  draws  a  box  around  the  portion  of  the  plot  to  be 
expanded,  and  then  pauses.  After  the  carriage  return  key  is  pressed,  the  boxed  por¬ 
tion  shown  in  figure  B4  is  then  replotted  as  shewn  in  figure  B5. 

Because  NPLOT  =  3  was  specified,  the  next  response  will  be  calculated  at 
mass  2.  As  shewn  in  table  BVI,  the  modal  accelerations,  acceleration  mobility,  and 
acceleration  response  for  mass  2  are  output  to  file  8.  The  program  loop  continues, 
as  shown  by  figure  B6,  table  BVII,  and  figure  B7,  until  the  NPLOT  parameter  is  equal 
to  a  positive  or  a  negative  5;  this  causes  the  program  to  terminate.  Table  BVIII 
shews  the  contents  of  file  8,  which  was  saved  after  finishing  the  interactive 
session. 
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TABLE  BI.-  MODAL  DATA  FOR  UNCONSTRAINED  SPRING  MASS 


Mode 

Frequency,  Hz 

Generalized  mass 

Mode  shape 

1 

0.0 

10.0 

1 

1.0 

1.0 

2 

7.749 

1 .179 

I 

.051649 

-.1  5761  5 

3 

15.442 

2.  514 

B 

615 

1.0 

-.234066 

TABLE  BII. -  INPUT  DATA  FILE  7  FOR  SPRING  MASS  SYSTEM 


SPRING  MASS  SYSTEM 
1  3  0  1  I  1  I  0.0  20.0  0.5 

NATURAL  FREQUENCY 
0.0  7.749  15.442 

GENERALIZED  MASS 
10.0  1.179  2.514 

DAMPING 

0.0 

FORCE  AT  MASS  3 

(1.0, 0.0) 

MASS  3  -  FORCE  COORDINATE 
1.0  -.157615  -.234066 

MASS  1  -  RESPONSE  COORDINATE 
1.0  1.0  -.3615 

MASS  2  -  RESPONSE  COORDINATE 
1.0  0.051649  1.0 


TABLE  Bill.-  CONTROL  PARAMETER  DESCRIPTION  FOR  SPRING  MASS  SYSTEM 


SPRING  HASS  SYSTEM 


CONTROL  PARAMETERS  SET 
NZ  -  1  (RIGID  »0D 


NZETA 

fCS 

NF 

ySTRT 

VSTOP 

DELU 


1  (RIGID  DODV  MODES* 

3  (NODES  TOTAL* 

•  (NORMAL  MODES* 
t  (VISCOUS  DAMPING* 

1  (DAMPING  IS  CONSTANT* 

1  (FORCE  VECTOR  SETS* 

1  (FORCE  COORDINATES* 

.9  (STARTING  FREQUENCY  (HZ)  ( 
39.9  (STOPPING  FREQUENCY  (HZ)  ( 
.599  (FREQUENCY  STEP* 


•LANK  COMMON  STORAGE  ZZZ  REQUIRES  379  LOCATIONS 
NOU  COMPUTING  FORCED  RESPONSE 
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TABLE  BIV. -  PRINT  MENU  FOR  MASS  1  RESPONSE 


SPRING  MASS  SYSTEM 

APPLIED  FORCE  UECTORi  FORCE  AT  MASS  3 
RESPONSE  LOCATION!  MASS  1  -  RESPONSE  COORDINATE 

FOR  THE  NEXT  3  QUESTIONS  ENTER  1  OF  THE  FOLLOWING  PRINT  OPTIONS: 

OPTIONS 

OUTPUT  TO  SCREEN . 1 

OUTPUT  TO  TAPE8 . 2 

OUTPUT  TO  SCREEN  fc  TAPE 8 . 3  • 

NO  OUTPUT  AT  ALL . 4 

MODaL  ACCELERATIONS’  (DEFAULT  3) 

’  2 

ACCELERATION  MOBILITIES’  (DEFAULT  3) 

’  2 

RESPONSE  TABLE’  (DEFAULT  3) 

i  2 

ENTER  NUMBER  OF  LINES  PER  PAGE  FOR  TAPE8  (DEFAULT  40) 

’50  • 

ENTER  PRINT  STE°  (DEFAULT  1) 

’  1 


TABLE  BV.-  NAMELIST  VPLOT  INPUT  FROM  TERMINAL  KEYBOARD  FOR  MASS  1  RESPONSE  PLOT 


ENTER  SUPLOT  NAHELISTt 
?  SUPLOT  QREAL-l,NPL0T-3  • 


-SUPLOT 

0QREAL 

1. 

0QIMAG 

9. 

0AP1PL 

9 . 

0PHASE 

9, 

0NPLOT 

3, 

0USCALE 

9, 

0XMIN 

0.0. 

0XMAX 

.  1E+02 

0YMIN 

0.0, 

0YMAX 

.  1E+02 

0JLINE 

1, 

eisvn 

0, 

0ISIZE 

1, 

0LTYPE 

9. 

0XDIU 

.  IE +02 

evDiu 

CSEND 

.  IE+02 

•  4 


I 


1  7 


NEED  GRIDS’  (V/N) 
?  N 
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TABLE  BVI . -  PRINT  MENU  FOR  MASS  2  RESPONSE 


SPRING  NASS  SYSTEM 

APPLIED  FORCE  UECTORt  FORCE  AT  NASS  3 
RESPONSE  LOCATION!  NASS  2  -  RESPONSE  COORDINATE 

FOR  THE  NEXT  3  QUESTIONS  ENTER  1  OF  THE  FOLLOWING  PRINT  OPTIONS* 

OPTIONS 


OUTPUT  TO  SCREEN . 1 

OUTPUT  TO  TAPE8 . 2 

OUTPUT  TO  SCREEN  &  TAPE8 . 3 

NO  OUTPUT  AT  ALL . 4 


NODAL  ACCELERATIONS?  (DEFAULT  3) 

?  2 

ACCELERATION  MOBILITIES?  (DEFAULT  3) 

7  2 

RESPONSE  TABLE?  (DEFAULT  3) 

?  2 

ENTER  NUMBER  OF  LINES  PER  PAGE  FOR  TAPE8  (DEFAULT  40) 

7  SO 

ENTER  PRINT  STEP  (DEFAULT  1) 

?  1 


TABLE  BV II.-  NAMELIST  VPLOT  INPUT  MASS  2  RESPONSE  PLOT 


ENTER  8UPL0T  NANEUSTi 
?  8UPL0T  QREAL‘1, ISYN-1 ,NPL0T«5  • 


-8UPL0T 

OQREAL 

00 I MAG 

OAMPL 

8PHASE 

ONPLOT 

OUSCALE 

exniN 

OXNAX 

0VNIN 

OYMAX 

OJLINE 

OISYM 

oisize 

OLTYPE 

OXDIU 

OYDIU 

08END 


1. 

e. 

o« 

e, 

5. 

o, 

0.0, 

.  1E+02, 

0.0, 

.  1E+02, 

1, 

1, 

1, 

0, 

. 1E+02, 
. iE+02, 


NEED  GRIDS7  (Y/N) 
?  N 
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TABLE  BVIII.-  FILE  8  SAMPLE  OUTPUT  OF  SPRING  MASS  SYSTEM  RESPONSES 


**********  V  I  B  R  A  ********** 
INPUT  DATA  CARD  IMAGE  FILE 


CARD 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 


SPRING  MASS  SYSTEM 
1301111  0.0  20.0  0.5 

NATURAL  FREQUENCY 
0.0  7.749  15.442 

GENERALIZED  MASS 
10.0  1.179  2.514 

DAMPING 

0.0 

FORCE  AT  MASS  3 

(1.0, 0.0) 

MASS  3  -  FORCE  COORDINATE 
1.0  -.157615  -.234066 

MASS  1  -  RESPONSE  COORDINATE 
1.0  1.0  -.3615 

MASS  2  -  RESPONSE  COORDINATE 
1.0  0.051649  1.0 


•  ( 


SPRING  MASS  SYSTEM 


APPLIED  FORCE  VECTOR:  FORCE  AT  MASS  3 
RESPONSE  LOCATION:  MASS  1  -  RESPONSE  COORDINATE 

MODAL  ACCELERATION  MATRIX  A(NF,NM) 

. 1000E+00  -. 1337E+00  .3366E-01 

SPRING  MASS  SYSTEM 


APPLIED  FORCE  VECTOR:  FORCE  AT  MASS  3 
RESPONSE  LOCATION:  MASS  1  -  RESPONSE  COORDINATE 


ACCELERATION  MOBILITY  Y (NR, NF) 


N 

FREQ 

YR 

YI 

1 

0.000 

.  1000E+00 

0. 

2 

.500 

.  1005E  +  00 

0. 

3 

1.000 

. 1021E+00 

0. 

4 

1.500 

.1049E+00 

0. 

5 

2.000 

.  1090E+00 

0. 

6 

2.500 

.  1146E  +  00 

0. 

7 

3.000 

. 1222E+00 

0. 

8 

3.500 

. 1324E+00 

0. 

9 

4.000 

.  1461E  +  00 

0. 

10 

4.500 

. 1649E+00 

0. 

11 

5.000 

.  1914E  +  00 

0. 

12 

5.500 

. 2308E+00 

0. 

13 

6.000 

.  294  2E+00 

0. 

14 

6.500 

.4101E+00 

0. 

15 

7.000 

. 6  84  3E+00 

0. 

16 

7.500 

. 2070E+01 

0. 

17 

8.000 

-. 2077E+01 

0. 

18 

8.500 

- .  706 1E  +  00 

0. 

19 

9.000 

-.  4341E+00 

0. 

20 

9.500 

-. 3200E+00 

0. 

21 

10.000 

-. 2589E+00 

0. 

22 

10.500 

-.2225E+00 

0. 

23 

11.000 

-. 2001E+00 

0. 

24 

11.500 

-. 1868E+00 

0. 

25 

12.000 

-.  1806E  +  00 

0. 

26 

12.500 

-. 1811E+00 

0. 

27 

13.000 

-. 1 893E+00 

0. 

28 

13.500 

-. 2085E+00 

0. 

29 

14.000 

- .  24  81E  +  00 

0. 

30 

14.500 

-. 3380E+00 

0. 

31 

15.000 

-. 6452E+00 

0. 

32 

15.500 

.  44  2  8E  +  01 

0. 

33 

16.000 

.  4165E  +  00 

0. 

34 

16.500 

.1996E+00 

0. 

35 

17.000 

.  1237E+00 

0. 

36 

17.500 

. 8575E-01 

0. 

37 

18.000 

.  6  338E-01 

0. 

38 

18.500 

. 4885E-01 

0. 

39 

19.000 

. 3879E-01 

0. 

40 

19.500 

. 3150E-01 

0. 

41 

20.000 

. 2604E-01 

0. 
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TABLE  BV I I I . -  Continued 


SPRING  MASS  SYSTEM 

APPLIED  FORCE  VECTOR:  FORCE  AT  MASS  3 
RESPONSE  LOCATION:  MASS  2  -  RESPONSE  COORDINATE 

MODAL  ACCELERATION  MATRIX  A(NF, NM) 

. 1000E+00  6905E-02  -.9311E-01 

SPRING  MASS  SYSTEM 

APPLIED  FORCE  VECTOR:  FORCE  AT  MASS  3 
RESPONSE  LOCATION:  MASS  2  -  RESPONSE  COORDINATE 

ACCELERATION  MOBILITY  Y(NR,NF) 


N 

FREQ 

YR 

i 

1 

0.000 

.  1000E+00 

0. 

2 

.500 

.  1001E+00 

0. 

3 

1.000 

.  1005E+00 

0. 

4 

1.500 

.  1012E+00 

0. 

5 

2.000 

.  1021E+00 

0. 

6 

2.500 

.  1033E+00 

0. 

7 

3.000 

. 104  9E+00 

0. 

8 

3.500 

.  1068E+00 

0. 

9 

4.000 

. 10  92E+00 

0. 

10 

4.500 

. 1122E+00 

0. 

11 

5.000 

. 115  8E+00 

0. 

12 

5.500 

. 1205E+00 

0. 

13 

6.000 

. 1269E+00 

0. 

14 

6.500 

. 1 364E+00 

0. 

15 

7.000 

. 1547E+00 

0. 

16 

7.500 

. 2310E+00 

0. 

17 

8.000 

. 2237E-01 

0. 

18 

8.500 

. 995  9E-0 1 

0. 

19 

9.000 

. 121 2E+00 

0. 

20 

9.500 

. 1361E+00 

0. 

21 

10.000 

. 1500E+00 

0. 

22 

10.500 

. 1649E+00 

0. 

23 

11.000 

. 1 822E+00 

0. 

24 

11.500 

. 2033E+00 

0. 

25 

12.000 

. 2301E+00 

0. 

26 

12.500 

. 2658E+00 

0. 

27 

13.000 

. 3158E+00 

0. 

28 

13.500 

. 3916E+00 

0. 

29 

14.000 

. 5199E+00 

0. 

30 

14.500 

. 7  844E+00 

0. 

31 

15.000 

.  1647E+01 

0. 

32 

15.500 

-.  1237E+02 

0. 

33 

16.000 

-.  1268E+01 

0. 

34 

16.500 

-. 65  89E+00 

0. 

35 

17.000 

-. 4411E+00 

0. 

36 

17.500 

-. 3292E+00 

0. 

37 

18.000 

-. 2611E+00 

0. 

38 

18.500 

-.  2154E+00 

0. 

39 

19.000 

-. 1 826E+00 

0. 

40 

19.500 

-.  157  9E+00 

0. 

41 

20.000 

-. 1387E+00 

0. 
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COMPUTER  REQUIREMENTS 

The  program  VIBRA  (vibration  response  analysis)  is  written  in  FORTRAN  V  for  the 
Control  Data  Corporation  6600  series  of  computers  and  implemented  under  the  NASA/LaRC 
Network  Operating  System. 


Memory  Allocation  and  Auxiliary  Storage  Files 

VIBRA  uses  dynamic  storage  allocation  for  memory  management  of  various  size 
problems.  All  the  arrays  are  packed  into  an  array  in  blank  common  called  ZZZ  in  the 
main  program.  The  length  of  the  blank  common  is  set  by  a  statement  within  the  pro¬ 
gram.  The  amount  of  blank  common  memory  required  for  a  particular  problem  is  deter¬ 
mined  by  the  total  number  of  frequency  points,  number  of  natural  modes,  force  vector 
sets,  and  force  coordinate  groups.  The  last  section  in  appendix  A  contains  informa¬ 
tion  for  estimating  the  length  of  the  blank  common  ZZZ  array.  If  there  is  insuf¬ 
ficient  blank  common  memory  for  a  problem,  VIBRA  terminates  with  an  error  message 
stating  the  present  length  of  array  ZZZ,  the  problem  size,  and  the  additional  amount 
of  memory  required. 

In  addition  to  using  blank  common,  several  auxiliary  files  are  used  during  pro¬ 
gram  execution  for  temporary  storage  and  recall  of  intermediate  calculations.  Three 
temporary  storage  files,  11,  12,  and  13,  are  suitable  for  saving  after  program  execu¬ 
tion,  in  addition  to  output  file  8.  Table  Cl  gives  a  description  of  the  program 
auxiliary  files.  The  use  of  files  5,  6,  7,  and  8  is  demonstrated  in  appendix  B,  with 
input  requirements  for  file  5  and  file  7  described  in  appendix  A.  As  listed  in 
table  Cl,  files  11,  12,  and  13  contain  modal  acceleration,  acceleration  mobility,  and 
acceleration  response  matrices,  respectively.  These  files  also  contain  index  infor¬ 
mation  that  will  enable  the  user  to  read  them  directly  into  another  computer  program. 
The  format  and  a  sample  content  of  these  files  are  given  in  tables  CII  to  CVII. 


Graphics  Subroutines 

VIBRA  uses  37  subroutines  from  the  PLOT-10  graphics  package,  which  consists  of 
the  Terminal  Control  System  (TCS),  described  in  reference  3,  and  the  Advanced 
Graphing  II  (AG-II),  described  in  reference  4.  Because  of  the  LaRC  Network  Operating 
System  (NOS)  configuration,  several  TCS  subroutines  have  been  modified.  The  purpose 
of  the  modifications  is  to  make  the  subroutines  run  more  efficiently  on  the  LaRC 
NOS.  However,  substitutions  can  be  made,  using  equivalent  TCS  subroutines,  to  pro¬ 
duce  the  same  functions  as  the  LaRC  modified  subroutines.  The  LaRC  NOS  modified 
PLOT-10  routines  (ref.  5)  deal  with  FORTRAN  input/output  while  in  graphics  mode  and 
with  the  file  output  buffer.  Also,  one  AG-II  subroutine,  called  LINE,  was  renamed 
ZLINE  to  avoid  a  conflict  with  one  of  the  LaRC  graphics  system  subroutines.  Another 
LaRC  NOS  graphics  modification  involves  the  placement  of  a  common  block,  named  JTB, 
in  the  main  program.  This  common  block  passes  parameters  to  the  NOS  PLOT-10  system 
which  controls  the  PLOT-10  vector  file  and  the  output  buffer  file  before  PLOT-1 0 
input/ output.  The  way  these  items  are  handled  may  be  different,  depending  upon  the 
version  of  PLOT-10  subroutines  available  on  various  computer  systems. 


•  4 


•  4 


•  4 


•  4 
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Tables  CVIII  to  CX  provide  lists  of  all  the  graphic  subroutines  used  and  their 
locations  within  VIBRA.  The  last  section  gives  a  detailed  description  of  the  LaRC 
modified  PLOT-10  subroutines  and  the  TCS  subprograms  called  by  them. 


LaRC  NOS  PLOT-10  Subroutines 

FORTRAN  subroutines  have  been  written  for  the  LaRC  NOS  system  and  use  a  few 
subroutines  from  the  TCS  package.  (See  ref.  5.)  The  purpose,  usage,  method,  and 
subprograms  used  by  these  NOS  subroutines  are  described  in  this  section.  Also,  a 
brief  description  of  the  JTB  common  block  used  in  the  main  program  is  given. 


Subroutine  EPAUSE 


Purpose:  To  execute  a  frame  advance 


Usage:  CALL  EPAUSE 


Method:  Terminate  a  plot  frame  and  go  to  the  next  frame, 

the  following: 


A  call  to  EPAUSE  causes 


1.  Flushes  the  plot  vector  file  buffer(s) 

2.  Rings  the  terminal  bell 

3.  Pauses 

4.  Continues  execution  when  user  presses  return  key 

5.  Does  an  immediate  screen  erasure 

Subprograms  used:  TSEND,  ERASE,  BELL 


Subroutine  TPAUSE 

Purpose:  To  dump  the  plot  vector  file  buffer  and  pause  for  a  user  prompt 

Usage:  CALL  TPAUSE 

Me thod:  Flush  the  plot  vector  file  buffer  to  the  terminal  and  pause.  A  call  to 

TPAUSE  causes  the  following: 


1.  Flushes  the  plot  vector  file  buffer(s) 

2.  Rings  the  terminal  bell 

3.  Pauses 

4.  Continues  execution  when  user  presses  return  key 
Subprograms  used:  TSEND,  BELL 
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Subroutine  NOTX 

Purpose:  To  display  horizontal  alphanumeric  labels 

Usage:  CALL  NOTX  (IX,  IY,  LENCHR,  ISTR),  where 

IX,  IY  Starting  position  of  label  in  screen  coordinates 

LENCHR  Number  of  characters  in  label  (maximum  of  100  characters) 

ISTR  Character  array  name  in  which  string  is  stored 

Method :  The  PLOT-10  label  routines  (HLABEL,  VLABEL,  AND  NOTAT)  require  labels  to  be 

stored  as  one  character  per  word,  ASCII  decimal  code.  This  is  an  inefficient  way 
to  store  labels  under  NOS.  NOTX  allows  user  to  store  labels  as  10  characters  per 
word,  left-justified  display  code.  The  labels  are  converted  from  display  code  to 
ASCII  and  are  displayed  by  the  PLOT-10  NOTAT  routine. 

Subprograms  used:  KAM2AS,  NOTAT 


Subroutine  NOTY  • 

Purpose:  To  display  vertical  alphanumeric  labels 

Usage :  CALL  NOTY  (IX,  IY,  LENCHR,  ISTR),  where 

IX,  IY  Starting  position  of  label  in  screen  coordinates  • 

LENCHR  Number  of  characters  in  label  (maximum  of  100) 

ISTR  Character  array  name  in  which  string  is  stored 

Method :  Same  as  NOTX  0 

Subprograms  used:  KAM2AS,  NOTAT 


Subroutine  DRASE 

Purpose:  To  erase  FORTRAN  output  from  the  Tektronix  screen 

Usage:  CALL  DRASE 

Method:  DRASE  should  be  used  (instead  of  ERASE  OR  NEWPAG)  on  NOS  to  erase  FORTRAN 

output  from  the  screen.  This  subroutine  does  not  write  an  erase  to  the  mass 
storage  plot  vector  file. 

Subprograms  used:  None 


•  1 
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JTB  Common  Block 

Ihe  JTB  common  block  is  located  in  the  main  program  only.  Default  values  have 
been  taken  from  all  the  common  block  parameters  except  for  JREQ. 

Purpose:  To  control  the  plot  vector  file  and  output  buffer 

Usage:  COMMON/ JTB/NFR,  JREQ,  IBAUD,  HDR,  I JO,  TFAC,  N FLUSH,  1JTB(3) 

NFR  =  0  (default) 

Pause  after  a  call  to  NFRAM  and  EPAUSE. 

Execution  continues  when  user  presses  return  key. 


JREQ 


IBAUD 


=  2 

PLOT-10  vector  file  is  assigned  to  terminal  and  also  written 
to  file  SAVPVF  for  saving. 

=  120  (default) 

Terminal  baud  rate  in  characters/second 


HDR  Used  internally 

I JO  Used  internally 


TFAC  =  (4014  terminal  length) /I  024 

N FLUSH  =  0  (default) 

Flush  the  user's  output  file  buffer  before  PLOT-10  input/output 


IJTB( 3) 


Reserved  for  future  use 


APPENDIX  C 


TABLE  Cl.-  PROGRAM  AUXILIARY  FILES 


File 

Function 

5 

Input  from  keyboard 

6 

Output  to  terminal  screen 

Input  data  file 

Output  (results  given  in  tabular  form) 

Modal  acceleration  matrices 

Acceleration  mobility  matrices 

Acceleration  response  matrices 

»  4 


TABLE  Cl I.-  AUXILIARY  FILE  11  FORMAT 


Card  image 

(a) 

Variable 

Description 

Format 

1 

NP 

Number  of  problems 

(13) 

NR 

Number  of  responses 

2 

1 

NF 

Number  of  forces 

>  (314) 

NM 

Number  of  modes 

) 

3 

PT 

Problem  title  card 

( A80) 

4 

FVT 

Force  vector  title 

( A80) 

5 

RCT 

Response  coordinate  title 

(A80) 

6 

MAT 

Modal  acceleration  heading 

( A80 ) 

7 

A(NR,NF,NM) 

Modal  acceleration  matrix 

1  0(  1  X,  El  0.4) 

aIf  there  is  more  than  one  problem  being  solved,  FILE  11  will  contain 
card  images  2  to  7  sequentially  repeated  NP  times.  Therefore,  card  image  1, 
which  contains  the  NP  variable,  defines  the  outer  loop  index.  The  index 
variables  on  card  image  2  set  up  the  inner  loop  for  reading  NR  groups  of  modal 
acceleration  matrices  calculated  for  one  problem. 


>  i 
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TABLE  CIII.-  FILE  11  CONTENTS  FROM  SAMPLE  PROBLEM  IN  APPENDIX  B 


2  13 

SPRING  MASS  SYSTEM 

APPLIED  FORCE  VECTOR:  FORCE  AT  MASS  3 
RESPONSE  LOCATION:  MASS  1  -  RESPONSE  COORDINATE 
MODAL  ACCELERATION  MATRIX  A(NF,NM) 

.  1000E+00  1337E+00  .3366E-01 

SPRING  MASS  SYSTEM 

APPLIED  FORCE  VECTOR:  FORCE  AT  MASS  3 
RESPONSE  LOCATION:  MASS  2  -  RESPONSE  COORDINATE 
MODAL  ACCELERATION  MATRIX  A (NFr NM) 

. 1000E+00  6905E-02  -.9311E-01 
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TABLE  CIV.-  AUXILIARY  FILE  12  FORMAT 


Card  image 

(a) 

Variable 

Description 

Format 

1 

NP 

Number  of  problems 

(13) 

|  NR 

Number  of  responses 

) 

2 

<  NF 

Number  of  forces 

(  (314) 

(  NM 

Number  of  modes 

) 

3 

IFN 

Number  of  frequency  steps 

(13) 

4 

PT 

Problem  title  card 

( A80 ) 

5 

FVT 

Force  vector  title 

(A80) 

6 

RCT 

Response  coordinate  title 

( A80 ) 

7 

AMT1 

Acceleration  mobility  title 

(A80) 

8 

AMT2 

Mobility  heading 

( A80) 

l  NSTEP 

Frequency  step  number 

)(1X,  13,  IX,  F8.3, 

9 

<  W(NR, NSTEP) 

Forcing  frequency 

?  1  2(1  X,  El  0.  4)  ) 

(YDD(NR, NSTEP, NF) 

Acceleration  mobility  matrix 

aIf  there  is  more  than  one  problem  being  solved,  FILE  12  will  contain  card 
images  2  to  9  sequentially  repeated  NP  times.  Therefore,  card  image  1,  which  con¬ 
tains  the  NP  variable,  defines  the  outer  loop  index.  The  index  variables  on  card 
image  2  set  up  the  inner  loop  for  reading  NR  groups  of  modal  acceleration  mobilities 
calculated  for  one  problem. 
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TABLE  CV.- 


FILE  12  CONTENTS  FROM  SAMPLE  PROBLEM  IN  APPENDIX  B 


1 

2 

4! 


SPRING  MASS  SYSTEM 

APPLIED  FORCE  VECTOR :  FORCE  AT  MAS:"  3 
RESPONSE  LOCATION:  MASS  1  -  RESPONSE 


ACC  ELF.  RAT  I  '• 

N  MOB  I L  TTY 

Y '  NR , 

N 

FREQ 

YR 

YI 

1 

0.000 

. 1000E+00 

0. 

.  500 

. 1 C05E+C0 

0. 

3 

1.000 

.  1071E+C-C 

0. 

4 

1.500 

. 1049F+C0 

0. 

5 

2.000 

. 1 09CE+0C 

0 . 

6 

2.500 

. 1 146E+00 

0 . 

7 

3.000 

. 1 222E*00 

0. 

8 

3. 500 

. 1324E+00 

0. 

9 

4.000 

. 1 46 1 E  *00 

0. 

10 

4.500 

.  16491:  *00 

0. 

11 

5.000 

. 1 9] 4E  *00 

0. 

12 

5.500 

.  230  8E  +  C0 

0. 

n 

6.000 

. 2942E-0C 

0. 

14 

6.500 

. 4  I C 1 E*QC 

0. 

15 

7.000 

.6843E-C0 

0. 

16 

7  .  500 

.  2070E*0 1 

0. 

17 

8.000 

-.2077E-01 

0. 

13 

8.500 

-.  7061E*G0 

0. 

19 

9. COO 

-.  434  1  E*00 

0. 

20 

9.500 

32C0E*0n 

0. 

21 

10.000 

2589E+00 

0. 

22 

10.500 

-. 2225E+00 

0. 

21 

11.000 

2001E+00 

0. 

24 

11.500 

- .  1  86  8E  +  00 

0. 

25 

12.000 

-.  1806E  +  00 

0. 

26 

12.500 

1  81 1E  +  00 

0. 

27 

13.000 

I  893  E  +  C0 

0. 

28 

13.50C 

-. 2085E*00 

0. 

29 

14.000 

- .  24  81 £  +  00 

0. 

30 

14.500 

-.  3  3  80E+00 

c. 

31 

15.000 

- .  64  5  2E  +  00 

0. 

32 

15.500 

. 4428E+01 

c. 

33 

16.000 

.  41 65E+00 

0. 

34 

16.500 

.  1996E  +  00 

0. 

35 

17. COO 

.  1237E  +  00 

0. 

36 

17.500 

. 857  5 E-0 1 

0. 

37 

18.000 

. 6  3  3  8E-0 1 

0. 

38 

18.500 

.  4  885E-0 1 

0. 

39 

19.000 

. 3879E-01 

0. 

40 

19.500 

.  3150E-01 

c. 

41 

20.000 

.  2604E-01 

0. 

41 

THING  MA: 

'.S  SYSTEM 

APPLIED  FORCE  VECTOR:  FORCE  AT  MASS 
RESPONSE  LOCATION:  MASS  2  -  RESPONSE 


ACCELERATION 

MOBILITY 

Y(NP,NF 

N 

FREQ 

YR 

YI 

1 

0.000 

. 1000E*00 

c . 

2 

.500 

. 1 00 1 E+00 

0. 

3 

1.000 

. 1005E+00 

0. 

4 

1.500 

. 1012E*00 

0. 

5 

2.000 

. 1G  2 1 E  +  00 

0. 

6 

2.500 

. 1033E+00 

0. 

7 

3.000 

. 104  9E  +  00 

0. 

8 

3.500 

. 106  SE  +  00 

0. 

9 

4.000 

. 1092E+00 

0. 

1C 

4.500 

. 1 1 22E+00 

c. 

11 

5.000 

. 2158E*00 

0. 

12 

5.50C 

.  1  2  0  5  E  ♦  0  0 

0. 

13 

6.000 

.  1  2  6  9  E  •  0  0 

0. 

14 

6.500 

. 1 364F  *00 

0. 

15 

".000 

.154^**00 

0. 

16 

7 . 500 

.  2  3 1  0  E  *  0  0 

0. 

17 

8. COO 

.2237F-01 

0 . 

1  8 

8.500 

. 9959E-Q1 

0. 

19 

9.000 

. 1 21 rt • 00 

0  . 

20 

9.  e0C 

. 1 36 I E*00 

c. 

21 

:o  .oor 

.  1  5  0  0  F  -  0  0 

0. 

2  2 

10.  =  00 

.  1649 E *00 

0  . 

2  3 

l :  .o oo 

.  1  92  2F.-C0 

0. 

.4 

1  1  .  500 

. 20 i If *00 

0. 

12.000 

. 2  30  IE *00 

G. 

26 

:  2  .  500 

.  26  5  8F.  *  CO 

0  . 

2" 

1  3.000 

.  J15  8F.*C0 

0 , 

;•  8 

: 1 . 5oo 

.  3  ^  1  6  E  *  0  0 

0  . 

2  o 

.  4  .  n  r.o 

. r  1  9  Q  p  *  0  0 

n . 

».;• 

'4.5  OH 

.  ’S44L*C0 

o  . 

i r .  r.  o  o 

.  164’F.*0l 

c. 

<  2 

: . 5 on 

.  I  2 J7F*P2 

f: . 

»  7 

:  f .  oor 

. :  .'6  pi-c ] 

■4 

;  f  . or. 

.  *i  5  9  O  F  *  0  0 

0  . 

:  7  .  t :• 

.44’.  IF  *00 

<(, 

o 

.  *r.c; 

c . 

3 ' 

\  9 .  :-(i> 

. 2* I ! F *00 

r, . 

*  8 

1  a .  =  o  r. 

. 2154E+00 

fl . 

<5 

:  Q .  r  f: 

.  I  82 6 F »CO 

o. 

40 

1  <1 .  •" 

.  ’.  6 -Of.  GO 

0  . 

4  1 

.  1  1  3  7  E  *  0  0 

0  . 

COORDINATE 


COORDINATE 


Reproduced  from 

bed  available  copy. 


•  < 


•  < 


•  < 


•  « 
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TABLE  CVI.-  AUXILIARY  FILE  13  FORMAT 


Card  imaqe 

Variable 

Description 

Format 

(a) 

1 

NP 

Number  of  problems 

(13) 

2 

NR 

Number  of  responses 

(13) 

3 

IFN 

Number  of  frequency  steps 

(13) 

4 

PT 

Problem  title 

( A80 ) 

5 

FVT 

Force  vector  title 

(A80) 

6 

RCT 

Response  coordinate  title 

( A80) 

7 

ART 

Acceleration  response  heading 

(A80) 

8 

(  NSTEP 

<  W(NR, NSTEP) 

( QDD( NR, NSTEP) 

Frequency  step  number 

Forcing  frequency 

Complex  acceleration  response 

)( 1  X,  13,  IX,  F8.3, 

I  2( 1 X,  El  0.  4)  ) 

aIf  there  is  more  than  one  problem  being  solved,  FILE  13  will  contain  card 
images  2  to  8  sequentially  repeated  NP  times.  Therefore,  card  image  1,  which  con¬ 
tains  the  NP  variable,  defines  the  outer  loop  index.  The  index  variables  on  card 
image  2  set  up  the  inner  loop  for  reading  NR  groups  of  acceleration  response  table 
calculated  for  one  problem. 


< 
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TABLE  CVII.- 


FILE  1 3  CONTENTS  FROM  SAMPLE  PROBLEM  IN  APPENDIX  B 


ft 


41 

SPRING  MASS  SYSTEM 

APPLIED  FORCE  VECTOR:  FORCE  AT  MASS  3 
RESPONSE  LOCATION:  MASS  1  -  RKSPONSI  COORD 1 NAT F 


N 

FREC 

OR 

V 

1 

n.ooo 

. 1C00E+C0 

0. 

2 

.  500 

. 1005 E+0C 

0. 

3 

1 .000 

. 1021E+00 

0. 

4 

1 . 500 

. ! 04  9E  *00 

0. 

5 

2.000 

. 1 0  9CE*00 

ft. 

6 

2.500 

. 1 146E*00 

0. 

7 

3.000 

. 1222E-OC 

c. 

a 

3.500 

. 1324E*00 

0. 

9 

4.000 

. 1461E+00 

0. 

10 

4.500 

. 1649E*00 

0. 

n 

5.000 

. 1 91 4E  *00 

0. 

12 

6.500 

. 2308E+00 

0. 

1  3 

6.000 

.  294  2E*00 

0. 

14 

6.500 

. 4101 E+C0 

0. 

IS 

7.000 

. 6843E+00 

0. 

16 

7.600 

. 2070E-C1 

c. 

17 

8.000 

2077 E * 0 1 

0. 

18 

3.500 

-.7061 E  +  00 

0. 

19 

9.000 

4341E+00 

0. 

20 

9.500 

3200E  +CC 

0. 

21 

10.000 

-. 2589E*00 

0. 

22 

10.500 

-.  2225E*00 

0. 

23 

11.000 

2001 E+00 

0. 

24 

11.500 

-.  1  868E  +  00 

o. 

25 

12.000 

-.  1  80  6  E ♦ 00 

0. 

26 

12.500 

-.1 81 1 E+00 

0. 

27 

13.000 

-. 1893E+00 

0. 

28 

13.500 

-. 2085E+00 

0. 

29 

14.000 

-.  24 81  E  +  00 

0. 

30 

14.500 

-. 3380E+00 

0. 

31 

15.000 

- . 6452E+00 

0. 

32 

15.500 

. 44  2  8E  +  0  1 

0. 

33 

16.000 

. 4 1 65E+00 

0. 

34 

16.500 

. 1996E+00 

0. 

?5 

17.000 

.  1  237E*00 

0. 

36 

17.500 

.  857  6E-01 

0. 

37 

18.000 

.  6  3  3  8E-0 1 

0. 

38 

18.500 

. 4  886 E -0 1 

0. 

39 

19.000 

. 3879E-01 

0. 

40 

19.500 

. 3150E-01 

0. 

41 

20.000 

. 2604E-01 

0. 

SPRING  MASS  SYSTEM 

APPLIED  FORCE  VECTOR:  FORCE  AT  MASS  3 
RESPONSE  LOCATION:  MASS  2  -  RESPONSE  COORDINATE 


FREQ 

OP 

01 

1 

0.000 

.  1000E  +  00 

0. 

2 

.  600 

.  1001 E  +  00 

0. 

2 

1.000 

.  1005E  +  00 

0. 

4 

1.600 

.  101 2E  +  00 

0. 

5 

2.000 

.  10  21 E+00 

ft. 

6 

2.500 

.  1O13F  +  C0 

0. 

7 

3.000 

.  1 C  4  9  E  ♦  0  0 

0. 

? 

3.500 

. 1G68E+00 

0. 

9 

4 .000 

.  1 0  92E  *00 

0. 

1  o 

4 . 5  0  0 

. 1 1 22E+00 

n . 

11 

5.000 

. 1 16  HE -00 

0. 

!  ?. 

6.500 

. 2206E*00 

0. 

1  3 

6.000 

.  1  2  6  9  E  *  0  0 

0. 

14 

6.500 

.  1  3  6  4  E  *  0  0 

c. 

16 

7.000 

.  1  54  - E ♦ 00 

0. 

16 

-.600 

.231  OF *00 

0  . 

P 

ft.  COO 

. 2237 K -01 

0. 

1  R 

P.  600 

. 9969E-0I 

0. 

I  9 

9.0  00 

. 1 21 2E+00 

0. 

70 

9. 600 

. 1 361 E*C0 

0. 

10.000 

. 1 500E ♦ 00 

0. 

10.500 

. 1649E*C0 

0. 

1  1.000 

.  1  ■* 2 7 E  -00 

0. 

.4 

1  1  .  500 

.  20  3  3E *00 

0. 

26 

1 2.000 

.  7  70  IK  +  00 

ft. 

76 

I  2  .  :  ft.  ft 

. 26  5  8E*00 

0. 

->  ■> 

1 3.(1  no 

. 3 1 5  ftp +00 

ft . 

'  * 

1 3.500 

.  3  9  3  6  E  +  0  0 

0  . 

,  9 

■4 .000 

.  6’  9  9E1  +00 

ft. 

*r 

.  4  on 

.  7  ,R4  4  K  +  on 

0  . 

• ) 

1 r  .  ■ :  0  o 

. 164’E+Ol 

0. 

1 r. .  6 or. 

-.  1  2  77E  +  02 

0. 

I  i 

1 6 .  i.or 

- . 1 26  PE +  01 

0. 

'4 

1  +  . 1 0r 

-  .  6  6  «  9  F  ♦  0  0 

0 . 

C 

1  • .  n  o  r. 

- . 441 i r +nn 

0. 

1  (-. 

1  • . !-  o  r 

1  292F  +  00 

0 . 

1  ft .  o  ft 

-.261  IE +00 

ft . 

1  P 

I  6.  50  0 

7 1  5=  4  1  +00 

0 . 

9 

.  9 .  f:pr- 

. 1 S26E+C0 

ft . 

4r 

: q .  r,.o 

;  qj  *r  n 

(!  . 

4 : 

; r .  ri.-f 

. J  1P-E*00 

0. 

» 
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TABLE  CVIII.-  PLOT-10  GRAPHICS  SUBROUTINES  USED  IN  VIBRA 


TCS  subroutines 

AG-II  subroutines 

ANMODE 

BINITT 

CHRSIZ 

CHECK 

DRAWA 

DLIMX 

ERASE 

DLIMY 

FINITT 

DSPLAY 

INITT 

FRAME 

MOVABS 

NPTS 

MOVEA 

SIZES 

TERM 

SLIMX 

TSEND 

SLIMY 

VCURSR 

STEPS 

SYMBL 

LaRC  subroutines 

XFRM 

XMFRM 

EPAUSE 

XMTCS 

TPAUSE 

XNEAT 

NOTX 

YFRM 

NOTY 

YMFRM 

DRASE 

YMTCS 

YNEAT 

ZLINE  (ZLINE  is  actually  LINE  in  AG-II) 

APPENDIX  C 


TABLE  CIX. -  PLOT- 10  SUBROUTINE  LOCATION  WITHIN  VIBRA 


PLOT-10  subroutine 

Location  in  VIBRA  subroutines 

ANMODE 

MAIN  GRAFIT  QUAD 4  ALABELS 

BINITT 

MAIN  TKSET  QUAD4  QRPV  QIPV  AMPV  PHPV 

CHECK 

GRAFIT  QRPV  QIPV  AMPV  PHPV  ZOOM 

CHRSIZ 

MAIN  TKSET  QRPV  QIPV  AMPV  PHPV 

DLIMX 

GRAFIT  ZOOM 

DLIMY 

GRAFIT  ZOOM 

DRASE3 

MAIN  PAUSE 

DRAWA 

ZOOM 

DSPLAY 

GRAFIT  QRPV  QIPV  AMPV  PHPV  ZOOM 

EPAUSE3 

MAIN  GRAFIT  QUAD 4  ALABELS 

ERASE 

TKSET  QUAD 4 

FINITT 

MAIN 

FRAME 

TKSET 

IN  ITT 

MAIN 

MOVABS 

MAIN  GRAFIT  QUAD 4  ALABELS 

MOVEA 

ZOOM 

NOTXa 

QRPV  QIPV  AMPV  PHPV  ALABELS 

NOTY3 

QRPV  QIPV  AMPV  PHPV  ALABELS 

NPTS 

GRAFIT  QRPV  QIPV  AMPV  PHPV  ZOOM 

SIZES 

TKSET 

SLIMX 

TKSET  QRPV  QIPV  AMPV  PHPV 

SLIMY 

TKSET  QRPV  QIPV  AMPV  PHPV 

STEPS 

TKSET 

SYMBL 

TKSET  QUAD 4 

TERM 

MAIN  TKSET  QUAD4 

TPAUSE3 

PAUSE 

TSEND 

GRAFIT  QUAD4  ALABELS 

VCURSR 

ZOOM 

XFRM 

TKSET 

XMFRM 

TKSET 

XMTCS 

TKSET 

XNEAT 

TKSET  QUAD 4 

YFRM 

TKSET 

YMFRM 

TKSET 

YMTCS 

TKSET 

YN  EAT 

TKSET  QUAD 4 

ZLINE° 

TKSET  QUAD 4 

aLaRC  NOS  PLOT-1 0  subroutine. 

^NOS  modified  name;  AG-II  name  is  LINE. 
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TABLE  CX.-  VIBRA  SUBROUTINES  THAT  USE  PLOT- 10  SUBROUTINES 


VIBRA  subroutine 

PLOT-10  subroutines 

MAIN 

ANMODE  BINITT  CHRSIZ  DRASE3  EPAUSE3 

FINITT  INITT  MOVABS  TERM 

PAUSE 

TPAUSE  DRASE3 

GRAFIT 

ANMODE  CHECK  DLIMX  DLIMY  DSPLAY 

EPAUSE3  FRAME  MOVABS  NPTS  SLIMX 

SLIMY  TSEND 

TKSET 

BINITT  CHRSIZ  ERASE  SIZES  STEPS  SYMBL 

TERM  XFRM  XMFRM  XMTCS  XNEAT  YFRM 

YMFRM  YMTCS  YNEAT  ZLINEb 

QUAD4 

ANMODE  BINITT  CHRSIZ  ERASE  MOVABS 

SYMBL  TERM  TSEND  XNEAT  YNEAT  ZLINEb 

QRPV 

QIPV 

AMPV 

PHPV 

BINITT  CHECK  DSPLAY  NOTX3  NOTY3  NPTS 

SLIMX  SLIMY 

ZOOM 

CHECK  DLIMX  DLIMY  DRAWA  DSPLAY  EPAUSE3 

NPTS  FRAME  MOVEA  VCURSR  XNEAT  YNEAT 

ALABELS 

ANMODE  MOVEABS  NOTX3  NOTY3  TSEND 

aLaRC  NOS  PLOT-1 0  subroutine. 

bNOS  modified  name;  AG-II  name  is  LINE. 
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TABLE  rr 


INPUT  DATA  FILE  7  FOR  CASE  1A 


TEST  CASE  1A  FOR  VIBRA  CONSTRAINED  SYSTEM - 

0  5  0  I  1  I  I  0.0  50.0  I. 

NATURAL  FREQUENCY  (HZ) 

6.01  14.05  22.13  27.51  42.47 


GENERALIZED 

MASS 

4.4  8  1.84 

6.  87 

2.  87  2.24 

VISCOUS  DAMPING 

.025 

FORCE  AT  NODE  1 

(1.0, 0.0) 
NODE  1  MODE 

SHAPE 

FOR  FORCE  AT  NODE 

1.00  -1.00 

1.00 

-.503  .00362 

RESPONSE  AT 

NODE  1 

1.00  -1.00 

1.00 

-.503  .00362 

RESPONSE  AT 

NODE  5 

.302  .273 

.926 

.317  1.00 

TABLE  III.-  CONTROL  PARAMETER  DESCRIPTION  FOR  CASE  1A 


TEST  CASE  1A  FOR  VIBRA  - CONSTRAINED  SYSTEM - 


CONTROL  PARAMETERS  SET 

N2  •  •  (RIGID  BODY  NODES* 

NN  »5  * NODES  TOTAL* 

NNORN  •  9  (NORMAL  NODES* 

MDANP  •  1  *UISC0US  DAMPING* 

NZETA  •  1  (DAMPING  IS  CONSTANT* 

NFS  •  1  (FORCE  VECTOR  SETS* 

NF  -1  (FORCE  COORDINATES* 

WSTRT  •  .9  (STARTING  FREQUENCY  (HZ)  * 

UST0P  -  59.9  (STOPPING  FREQUENCY  (HZ)  * 

DELU  •  1.999  (FREQUENCY  STEP* 

BLANK  COMMON  STORAGE  Z2Z  REQUIRES  3S1  LOCATIONS 

NOW  COMPUTING  FORCED  RESPONSE 
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6  00  E 


TABLE  VIII.-  INPUT  DATA  FILE  7  FOR  CASE  IB 


TEST  CASE  IB  FOR  VIBRA  CONSTRAINED  SYSTEM 

05  0  1  0  2  3  20.0  20.0  0. 

NATURAL  FREQUENCY  (HZ) 

6.01  14.05  22.13  27.51  42.47 

GENERALIZED  MASS 
4.48  1.  84  6.  87  2.  87  2.24 

VISCOUS  DAMPING 
.025  .04  .018  .02  .06 

♦FORCE  VECTOR  SET  1* 

(100,-50)  (-50,-200)  (150,150) 

♦FORCE  VECTOR  SET  2* 

(50,-50)  (250,-100)  (-100,50) 

NODE  1  MODE  SHAPE  FOR  FORCE  AT  NODE  1 
1.00  -1.00  1.00  -.503  .00362 

NODE  2  MODE  SHAPE  FOR  FORCE  AT  NODE  2 


.857  -.221 

-.934 

1.00 

-.0221 

NODE  3  MODE  SHAPE  FOR  FORCE 

AT  NODE 

.694  .298 

-.547 

-.489 

.0829 

RESPONSE  AT 

NODE  1 

1.00  -1.00 

1.00 

-.503 

.00362 

RESPONSE  AT 

NODE  2 

.  857  -.  221 

-.934 

1.00 

-.0221 

RESPONSE  AT 

NODE  3 

.694  .298 

-.547 

-.489 

.0829 

RESPONSE  AT 

NODE  4 

.470  .373 

.944 

.223  -. 

536 

RESPONSE  AT 

NODE  5 

.302  .273 

.926 

.317  1. 

00 

TABLE  IX.-  CONTROL  PARAMETERS  DESCRIPTION  FOR  CASE  IB 


TEST  CASE  IS  FOR  VIIRA  —CONSTRAINED  SYSTEM - 


CONTROL  PARAMETERS  SET 

NZ  •  9  (RIGID  BODV  NOSES* 

NN  ■  5  (NODES  TOTAL* 

NNORN  •  9  (NORMAL  NODES* 

WDANP  •  1  (VISCOUS  DAMPING* 

NZETA  •  9  (DAMPING  VARIES* 

NFS  •  2  (FORCE  VECTOR  SETS* 

NF  •  3  (FORCE  COORDINATES* 

USTRT  •  29.9  (STARTING  FREQUENCV  (HZ!  * 

UST0P  ■  29.9  (STOPPING  FREQUENCV  (HZ)  * 

DCLU  •  .999  (FREQUENCY  STEP* 

BLANK  COMMON  STORAGE  ZZZ  REQUIRES  79  LOCATIONS 

NOW  COMPUTING  FORCED  RESPONSE 
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TABLE  XII.-  ACCELERATION  RESPONSE  AT  20  Hz  FOR  NODES  3,  4  AND 


APPLIED  FORCE  VECTOR:  *F0RCE  VECTOR  SET  2* 
RESPONSE  LOCATION:  RESPONSE  AT  NODE  5 


< 


TABLE  XIII.-  INPUT  DATA  FILE  7  FOR  CASE  2 


OH  5  8  TAILBOOM  VIBRATION  RESPONSE 


6  12  1  0  0  1 
NATURAL  FREQUENCY 
. 99269222E-03 
.  15454463E-02 
. 21336986E+03 
GENERALIZED  MASS 
1.0 

DAMPING 

0  0  0  0  0  0 


1  75.00  : 

(HZ) 

10557  244E-02 
1 85  85053E-02 
216  21 83  5E+0  3 


0.015  0.016 


00  1.00 

.  1193  82  88E-02 
. 7  9266240E+02 
.  2967  820  9E  +  03 


. 14050625E-02 
. 7  9623  835E+02 
.  36  830  964E+0  3 


FORCE  VECTOR  1.0 

(-1.0, 0.0) 

FORCE  NODE  193 
.  3024  8530E+01 
.  941514  86E-01 
-. 4  8201405E+01 
RESPONSE  NODE  109 
.  9937  2428E+00 
-. 86123459E-01 
-.  1097  8691E+01 


1.0  LBS  DOWN  AT  NODE 


016  0.013  0.015  0.016 

193 


DIRECTION 
-.  37  951904E+01 
-.  36754599E+01 
. 21011564E+01 
1  DIRECTION 
-.  26  37  237  8E+01 
-.  984397  80E+00 
. 4  995  36  24E+00 


.  15550316E+01 
5  26  81 4  21 E+01 
141 867  34E+01 

17  557  607E+01 
3  8250634E+01 
50664644E+00 


. 25379737E+01 
-.  26225012E+01 
-.  16998924E  +  00 

-.  12423869E+00 
. 19295327E+01 
. 24562881E+01 


TABLE  XIV.-  CONTROL  PARAMETER  DESCRIPTION  FOR  CASE  2 


0H68  TAILBOON  VIBRATION  RESPONSE 


CONTROL  PARANETERS  SET 


S  SRIGID  BODY  NODESS 
18  INODES  TOTALS 
1  SORTHONORNAL  NODESS 
0  SSTRUCTURAL  DANPINGS 
9  SDANPING  VARIES! 

1  SF0RCE  VECTOR  SETS! 

1  SFORCE  COORDINATES! 

75. 0  SSTARTING  FREQUENCY  (HZ)  S 
3OO.0  SSTOPPINC  FREQUENCY  (HZ>  S 
1.000  SFREQUENCY  STEPS 

BLANK  CONNON  STORAGE  ZZZ  REQUIRES  1443  LOCATIONS 
NOW  CONPUTING  FORCED  RESPONSE 
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Figure  1  .  - 


Flow  chart  of  VIBRA  program 


rtocHHr'oas  r-  x>  m 


APPLIED  FORCE  UECTORI  force  at  NODE  1 
RESPONSE  LOCATION*  RESPONSE  AT  NODE  X 


Figure  3.-  Acceleration  response  components  at  node  1  for  case  1A 


Figure  8.-  Expanded  plot  of  real  acceleration  component  from  figure  7 


APPLIED  FORCE  VECTOR* 
RESPONSE  LOCATION I 


FORCE  VECTOR  l.S  LIS  DOUN  AT  NOSE  183 
RESPONSE  NODE  1SB  1  DIRECTION 
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